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AbSttACt 

Me  consider  the  stability  to  thermal  perturbation  of  static 
models  of  coronal  loops.  The  effects  of  cool,  radiatively  stable 
material  at  the  loop  base  are  included.  The  linear  stability 
turns  out  to  be  sensitive  only  to  the  boundary  conditions  assumed 
on  tne  velocity  at  the  loop  base.  The  question  of  the 
appropriate  boundary  conditions  is  discussed,  and  we  conclude 
tnat  tne  free-surtace  condition  (the  pressure  perturbation 
vanishes),  rather  than  the  rigid-wall  (the  velocity  vanishes),  is 
relevant  to  the  solar  case.  Me  find  that  the  static  models  are 
thermally  unstable,  with  a  growth  time  of  the  order  of  the 
coiunai  cooling \time.  The  physical  implications  of  these  results 
for  the  solar  corona  and  transition-region  are  discussed.  " 
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The  linear  theory  for  the  theraal  stability  of  static 
coronal- loop  aodels  has  been  investigated  by  a  nuaber  of  authors, 
bur  wlcn  differing  conclusions.  In  his  original  work,  Antiochos 
(1979)  concluded  that  the  static  aodels  were  theraally  unstable. 
Siaiiar  results  were  obtained  by  Hood  and  Priest  (1980). 
However,  Chiuderi,  Einaudi,  and  Torricelli-Ciaaponi  (1981),  Craig 
ana  McClyaont  (1981),  and  McClyaont  and  Craig  (1981a, b,c)  have 
found  that  the  aodels  are  either  stable  or  that  the  growth  rates 
for  instability  are  too  saall  to  be  physically  significant.  The 
reason  for  this  difference  in  the  results  of  the  two  sets  of 
auenors  is  in  tneir  treatment  of  the  base  of  the  loop  aodels. 
Antiochos  (1979)  and  Hood  and  Priest  (1980)  have  not  included 
cool  aateriai,  T  -  105  K,  in  their  aodel,  for  which  the  fora  of 
radiative  loss  curve  (e.g.,  Rayaond,  Cox,  and  Saith  1978)  favors 
linear  stability  (Field  1965).  In  addition,  Antiochos  has  consi¬ 
dered  only  perturbations  with  a  vanishing  first-order  heat  flux 
at  the  base.  Chiuderi,  Einaudi,  and  Torricelli-Ciaaponi;  and 
Craig  and  McClyaont  argue  that  the  growth  rates  for  instability 
are  very  sensitive  to  these  assuaptions  so  that  the  aodels  can  be 
effectively  stabilised  by  either  including  soae  chroaospheric 
aateriai  or  by  changing  the  boundary  conditions  on  the  teapera- 
ture  perturbation;  in  particular,  the  aodels  can  be  aetastable 
(sero  growth  rate)  if  the  teaperature  perturbation,  instead  of 
the  heat-flux  perturbation,  is  assuaed  to  vanish  at  the  base. 


It  is  clear,  therefore,  that  the  proper  treatment  of  the 
base  is  critical  for  determining  to  what  extent  the  aodels  are 
eicuer  stable  or  unstable.  In  the  next  section  we  discuss  the 
question  of  the  boundary  conditions  in  detail. 

ii.  gfittmb&tlan  Equations  and  Boundary  Conditions 

Proa  the  fora  of  the  perturbation  equations,  e.g.,  Antiochos 
(1979),  Field  (1965),  it  is  evident  that  there  are  four 
independent  spatial  derivatives  in  the  problea  and,  hence,  only 
four  boundary  conditions  can  be  specified.  This  is  to  be 
expected  since  the  full  nonlinear  equations  also  require  four 
spatiax  boundary  conditions  for  a  unique  solution  (e.g., 
Richtayer  and  Norton  1967).  Bence,  at  each  end  of  the  loop  two 
conditions  aust  be  specified.  One  condition  generally  describes 
the  theraal  properties  of  the  base;  for  exaaple,  the  loop  base 
may  act  as  a  thermal  bath  so  that  the  temperature  perturbation 
vanishes  there,  Tlb  -  0;  or  it  aay  act  as  a  thermal  insulator  so 
tnat  the  perturbation  heat  flux  vanishes  there,  P^b  »  0.  The 
other  condition  describes  the  inertial  properties  of  the  base; 
for  exaaple,  it  aay  act  as  a  rigid  wall  so  that  the  velocity 
perturbation  vanishes,  v^b  ■  0,  or  a  free  surface,  so  that  the 
pressure  perturbation,  Plb  ■  0.  Depending  on  the  physical  situa¬ 
tion,  either  of  the  inertial  conditions  may  be  assumed,  but 
clearly  not  both.  Assuming  both  conditions  is  physically  incon¬ 
sistent;  it  implies  that  the  loop  base  is  both  rigid  and  free. 
Also,  it  is  mathematically  incorrect  since,  along  with  the 
temperature  condition,  it  is  equivalent  to  imposing  six  indepen¬ 
dent  boundary  conditions  on  a  problem  that  admits  only  four. 
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However#  previous  authors  (Antiochos  1979|  Chiuderi# 
Einaudi#  and  Torricelli-Ciamponi  1980}  and  Craig  and  McClymont 
1981;  have#  in  fact#  iaposed  both  inertial  boundary  conditions  on 
the  perturbation.  In  their  case#  they  were  able  to  find 
solutions  to  the  overconstrained  problem  because  they  considered 
only  the  restricted  class  of  models  in  which  gravity  is 
neglected;  and  the  loop  area  and  the  coronal  heating  do  not  have 
any  spatial  variation.  Under  these  siaplifications#  the  static 
equations  are  autonoaous#  and  the  equilibriua  models  are 
symmetric  about  the  loop  apex.  Hence#  the  solutions  to  the 
first-order  equations#  i.e.#  the  normal  modes  of  the  loop#  are 
either  purely  symmetric  or  antisymmetric  about  the  apex.  But# 
for  the  simplified  system  with  no  gravity#  the  only  possible 
antisymmetric  solution  to  the  force  equation  is  that  vanishes 
identically.  Hence#  all  the  antisymmetric  solutions#  irrespec¬ 
tive  of  their  boundary  conditions#  will  also  trivially  satisfy 
tne  free-surtace  conditions}  in  particular#  all  antisymmetric 
modes  for  the  rigid-wall  conditions  will  satisfy  the  overcon¬ 
strained  problem. 

Previous  authors  have  used  this  result  to  exclude  the 
symmetric  modes  as  valid  solutions.  The  argument  used  is  that 
tne  antisymmetric  modes  are  more  physical  because  all  the 
perturbed  quantities#  T^#  v^#  and  P^  vanish  at  the  boundary. 
But  it  is.  evident  from  the  discussion  above  that  this  argument  is 
spurious  since#  in  general#  no  mode  can  satisfy  all  these  condi¬ 
tions.  Por  a  realistic  equilibrium  model  that  has  no  special 
symmetry  properties#  the  modes  will  be  neither  symmetric  nor 
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antisymmetric.  Instead*  there  will  be  two  distinct  ssts  of  nodes 
corresponding  to  the  two  distinct  physical  conditions  of  either 
rigid  wells  or  free  surfaces*  end  no  node  will  satisfy  both. 
Bence*  within  each  set  ell  the  aodes  ere  equally  acceptable*  or 
unacceptable. 

The  question  remains  as  to  which  set  of  boundary  conditions 
best  represents  the  physical  situation  in  a  solar  loop.  Ideally* 
one  would  like  to  place  the  loop  boundary  sufficiently  deep  down 
in  the  ataosphere  that  conditions  there  do  not  affect  the 
stability  properties  of  the  corona  and  transition  region.  This 
can  be  done  for  the  theraal  boundary  conditions  by  placing  the 
loop  base  at  a  sufficiently  low  temperature.  Since  both  conduc¬ 
tion  and  radiation  decrease  rapidly  for  T  <  6  x  10*  K*  the 
tneraai  time  scale  at  the  base  can  be  made  to  be  very  long 
compared  to  the  time  scales  in  the  corona  or  transition  region 
simply  by  including  sufficiently  cool  material  at  the  base  of  the 
static  models.  We  expect  that  in  this  case  the  growth  rates  and 
tne  eigenfunctions  will  be  insensitive  to  the  thermal  boundary 
conditions. 

However*  it  is  not  possible  to  select  the  base  point  so  that 
the  growth  rates  and  eigenfunctions  are  independent  of  the  iner¬ 
tia*  boundary  conditions.  The  reason  for  this  is  that  the  sound 
travel  time  in  the  corona*  chromosphere  and  even  photosphere  is 
very  rapid  compared  to  typical  thermal  time  scales*  such  as  the 
coronal  cooling  time.  Bence*  we  expect  that  the  stability 
properties  of  the  static  models  will*  in  general*  be  sensitive  to 
the  inertial  boundary  conditions  assumed  for  the  perturbations; 
we  show  below  that  this  is  indeed  the  case. 
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We  art  left,  therefore,  with  having  to  dacida  which  is  tba 
propar  inartial  boundary  condition  at  tha  loop  basa.  We  baliava 
that  tha  corract  condition  is  tha  fraa  surf scat  Plb  ■  0,  and 
that  tha  rigid-wall  condition  is  not  applicabla  to  solar  loops. 
Assuming  tha  rigid-wall  condition  iaplias  that  at  soaa  laval  in 
tha  solar  ataosphara  a  trua  physical  boundary  exists  across  which 
no  aass  transfer  can  occur.  Ona  may  think  of  a  rigid  wall  as 
representing  aatarial  of  infinite  inertia  and,  hence,  imaobile. 
For  axaapla,  in  tha  loop  models,  the  magnetic  field  is  assumed  to 
provide  infinite  inertia  to  motions  perpendicular  to  the  field 
direction,  thereby  justifying  a  one-dimensional  loop  geometry. 
However,  in  our  model  the  plasma  itself  has  essentially  no 
inertia  since  that  is  the  assumption  implicit  in  dropping  the 
acceleration  terms.  Therefore,  by  assuaing  rigid-wall  conditions 
at  tha  base,  we  are  postulating  that  there  is  a  distinct  physical 
difference  between  tha  plasaa  in  the  loop  and  the  plasaa  below, 
i.e.,  the  inertial  properties  of  the  solar  ataosphere  must  change 
abruptly  at  some  level.  To  our  knowledge,  there  axists  no  such 
level. 

It  is  important  to  note  that  such  an  inertial  change  is  not 
equivalent  to  abrupt  density  or  pressure  changes,  which  do  occur 
in  the  atmosphere.  In  particular,  tha  chromosphere  may  be 
thought  of  as  a  narrow  interface  between  the  low  density  and 
pressure  corona,  and  tha  high  density  and  pressure  photosphere. 
Due  to  the  density  difference,  the  velocities  in  the  photosphere 
will  tend  to  be  very  small  compared  to  those  in  tha  corona.  For 
example,  tha  velocities  in  tha  photosphere  beneath  a  coronal  hole 


are  negligible  compared  to  solar  wind  speeds.  But  the  important 
quantity  is  the  mass  flux,  and  this  is  the  same  in  both  regions. 
Also,  due  to  the  pressure  ratio,  the  position  of  the  top  of  the 
photosphere  will  be  negligibly  affected  by  pressure  changes  in 
the  corona,  but  this  does  not  mean  that  the  top  of  the 
photosphere  acts  as  a  rigid  wall  with  respect  to  the  corona.  On 
tne  contrary,  material  moves  freely  from  the  photosphere  to  the 
corona  and  vice  versa.  If  the  top  of  the  photosphere  acted  as  a 
true  rigid  wall,  the  mass  flux  would  necessarily  vanish  there; 
and  there  would  be  no  steady-state  flows  such  as  the  solar  wind 
or  siphon  flow  into  sunspots. 

hi.  OLQMlh  SAtfiA  and  BigenfunctionA 

A.  Free-Surface  Nodes 

In  this  letter  we  calculate  the  growth  rates  and  eigen¬ 
functions  only  for  the  simplified  problem  that  has  been  discussed 
by  previous  authors;  more  realistic  models  will  be  considered 
in  a  subsequent  paper.  In  addition,  we  consider  only  the  so- 
called  thermally  isolated  static  model,  in  which  the  heat  flux  is 
assumed  to  vanish  at  the  loop  base  (Vesecky,  Antiochos  and 
Underwood  1979).  However,  we  investigate  in  detail  the  effects 
of  adding  cool  material  to  the  coronal  loop  by  including  a  model 
chromosphere,  and  we  include  in  the  analyses  the  first  symmetric 
mode,  which  we  have  omitted  previously.  For  comparison  we  calcu¬ 
late  the  growth  rates  for  both  the  free-surface  and  the  rigid- 
wall-boundary  conditions;  although  as  discussed  above,  we  believe 
tnat  the  free-surface  condition  is  the  physically  relevant  one. 


Under  the  assumptions  above,  and  neglecting  the  acceleration 
teraa  (Field  1965),  the  firat'-order  equations  can  be  written  in 


fora  aat 

C"  +  Q(x)  £  +  A  T0~3/2£-  P1  <.4XT0  +  G(x))/PQ  (1) 

where t 

€  -  T03/2  Tx  ,  (2) 

Q(x)  ■  (n0A  +  T0  6  A  T0  (e/n0  -  n0A  ) 

-  nQ  6/6ne  (e/no))/(10“6  n0T07/2)  ,  (3) 

G(x)  *  (e/n0  -  nQ  6/6  nc  (e  /no))/(10“6noTo)  ,  (4) 

A  -  -  107  k2  v/p0  .  .  (5) 


In  (1)  -  (5)  P  is  the  total  plasma  pressure;  n  is  the  electron 
density;  A  is  the  radiative  loss  coefficient  (e.g.,  Raymond,  Cox 
and  Saith  1976);  e  is  the  coronal  heating  function;  the 
coefficient  of  theraal  conduction  is  that  of  Spitser  (1962);  the 
subscript  *o"  refers  to  the  static  aodel;  k  is  Boltsaann's 
constant;  v  is  the  growth  rate,  i.e.  all  first-order  variables 
are  assuaed  to  vary  as  eVfc;  and  the  prime  indicates  differentia¬ 
tion  with  respect  to  the  independent  variable  x,  which 
corresponds  to  the  total  nuaber  of  electrons  initially  in  the 
loop  and  aeasured  froa  the  top  (see  Antiochos  1976  or  NcClyaont 
and  Craig  1981b)* 

Equation  (1)  aust  now  be  solved  for  the  various  boundary 
conditions*  The  theraal  boundary  conditions  are  straightforward 
to  iapleaentt  if  the  teaperaturi  perturb*  ion  i  assuaed  to 


vanish  at  the  base,  then  £b  -  0;  if  the  heat  flux  perturbation 
vanishes,  then  <S')b  ■  0.  This  is  the  main  advantage  to  using  a 
Lagrangian  formulation.  In  Bulerian  coordinates  the  boundary 
condition  of  a  vanishing  perturbation  heat  flux  has  an  integral 
form  in  general;  whereas  it  has  the  simple  differential  form 
above  in  Lagrangian  coordinates.  Note  that  there  is  a  physical 
difference  between  the  two  formulations.  In  the  Eulerian  case  we 
are  setting  the  boundary  at  a  particular  point  in  the  loop, 
wnereas  in  the  Lagrangian  case  we  are  setting  the  boundary  at  a 
particular  plasma  particle.  (Of  course,  both  formulations  are 
identical  in  the  case  of  rigid-wall  conditions).  Since  we  will 
show  below  that  the  stability  properties  are  insensitive  to  the 
particular  position  of  the  boundary,  as  long  as  it  is  somewhere 
in  the  chromosphere,  we  may  as  well  use  the  Lagrangian 
formulation  for  mathematical  convenience. 

The  free-surface  condition  is  also  straightforward.  Since 
Pj  is  a  constant  fc_  the  case  of  no  gravity,  the  free-surface 
condition  implies  that  the  righthand  side  of  (1)  vanishes. 
Hence,  we  are  left  with  a  standard  Sturm-Liouville  problem  as  in 
Antiochos  (1979),  which  we  solve  in  exactly  the  same  manner  as 
before.  The  rigid-wall  condition,  on  the  other  hand,  is  more 
difficult  to  implement.  The  simplest  procedure  is  to  express  it 
as  an  integral  constraint,  equivalent  to  Equation  (20)  in 
Antiochos  (1979).  Letting  the  position  of  the  base  points  be  ±X, 
(±L  in  Eulerian  coordinates),  we  find  that  the  perturbed 
pressure,  P. ,  can  be  determined  in  terms  of  the  eigenfunctions: 


dx 


(6) 


X 

px  -  (k/L)  j  i  T0"3/2 
-X 

Therefore,  for  the  rigid-wall  condition,  we  must  solve  the 
inhomogeneous  problem  (1),  subject  to  the  constraint  (6)  (cf. 
McClymont  and  Craig  1981c).  Note  that  this  is  not  a  standard 
Sturm-Liouville  problem;  hence,  none  of  the  standard  theorems 
(e.g.,  Morse  and  Feshbach  1953)  on  the  behavior  of  the 
eigenfunctions  and  eigenvalues  need  apply.  Since  for  this  prob¬ 
lem  tnere  are  two  parameters,  X  and  P^,  that  must  be  determined, 
we  do  not  use  a  "shooting”  technique  as  in  the  homogeneous  case 
to  obtain  solutions.  Instead,  we  consider  (1)  and  (6)  as 
defining  a  nonlinear  problem  for  X  andC/P^  and  use  basically 
Newton-Raphson  iteration.  Convergence  was  rapid,  -  20 
iterations,  for  all  the  cases  we  ran. 

There  are  two  functions  in  the  problem  that  must  be  speci¬ 
fied:  the  radiative  loss  coefficient,  A (T) ,  and  the  energy  input 
per  particle,  e  <n,T)/n.  For  the  former  we  use  a  smooth  analytic 
approximation  to  the  curve  derived  by  Raymond,  Cox  and  Smith 
(1976) .  The  stability  properties  of  the  models  are  not  sensitive 
to  the  exact  form  of  A(T)  as  long  as  this  form  is  such  that  it 
implies  radiative  instability  above  '  105  k  and  stability  below  ~ 
105  K. 

Since  the  mechanism  for  coronal  heating  is  not  known,  the 
form  of  the  energy  input  £  is  essentially  arbitrary.  In  our 
previous  work  we  found  that  the  stability  properties  are  also 
insensitive  to  the  exact  form  for  £  ,  at  least,  for  the  case 
where  this  form  is  a  power  law  dependence  on  n  and  T.  Bence,  we 
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simply  assume  that  in  the  cocona  and  transition  region,  i.e.,  T  - 
6  x  104K,  the  energy  input  per  particle  is  constant:  e  /n  • 
const.  However,  for  the  lower  temperatures  we  use  a  different 
form  for  the  energy  input  so  that  there  will  be  a  significant 
mass  of  material  at  "chromospheric"  temperature,  T  <  6  x  104  K. 
The  simplest  procedure  for  obtaining  a  thick  chromosphere  in  the 
static  models  is  to  adjust  the  energy  input  rate  at  low  tempera¬ 
tures  so  that  it  becomes  almost  equal  to  the  radiation  loss  rate 
(e.g.,  Craig,  Robb  and  Rollo  1981).  Therefore,  we  use  the 
following  form  for  the  energy  input: 

e/n  -  C0  tanh(x<T)>  +  nA  (1  -  tanh(x(T)»  (7) 

where 

X  (T)  -  Cx  <T/Tob  -  1  +  6>m  (8) 

and 

CQ,  Cj,  Tob,  m  and  5  are  constants.  Equations  (7)  and  (8)  imply 
that  for  5  <<  1:  e/n  ->  CQ  for  T  >  Tob  and  that  e/n  ->  nA  for  T  - 
Tob.  The  amount  of  material  at  Tob  becomes  arbitrarily  large  for 
arbitrarily  small  6 »  in  particular,  one  finds  a  chromospheric 
size  scale  *1/6  for  the  value,  m  »  3,  of  the  exponent  in 
Equation  (8).  Of  course,  this  is  a  highly  unrealistic  model  for 
the  true  solar  chromosphere.  However,  it  permits  us  to 
investigate  conveniently  the  effect  of  adding  cool,  radiatively 
stable  material  to  the  loop  base. 

ffe  have  analyzed  the  stability  properties  of  a  wide  range  of 
static  models.  In  Figure  1  we  plot  the  growth  rates  of  the  three 
lowest  free-surface  modes  for  a  series  of  static  models.  These 


Or  with  either 


modes  are  the  first  symmetric  modes#  ?i,£0p 
t3  base  *  0  or  Tl,base  “  0 •  and  the  first  antisymmetric  mode# 
Tl , top  “  °»  with  T^base  *  0.  The  static  models  all  have  almost 
identical  coronal  properties#  i.e.#  at  the  loop  apex  the  tempera¬ 
ture,  Tcor  ■  2.3  x  106  K#  the  density,  ncor  ■  5.5  x  109  cm"3#  and 
the  half-length#  L  -  1.7  x  109  cm.  The  models  were  generated  by 
fixing:  the  energy  input  rate  into  the  corona#  CQ  ■  10"32 
ergs/sec/particle  [Equation  (7)1;  the  chromospheric  temperature# 
Tob  “  3  x  104  K;  and  the  exponent  m  ■  3.  We  vary  either  the  base 
temperature  (Tb)  #  or  the  mass  of  the  base  (Tb  -  Tob  and  6 
varies) . 

The  growth  rates  in  Figure  1  are  expressed  in  terms  of  the 
conductive  cooling  time  of  the  coronal  plasma#  i.e.#  we  plot  X  > 
vtc#  where  ^  is  defined  as: 


T 


c 


5/2  PQ  L2 

10**5  m  7/2 
AU  Acor 


(9) 


For  the  particular  values  of  the  loop  parameters  above#  c  ■  1.4 
x  103  sec;  however#  is  insensitive  to  the  value  of  c.  Hence# 
Figure  1  should  be  valid  for  any  static  model  whose  parameters 
lie  within  the  observed  range  of  solar  values.  Note  that  for  all 
the  free-surface  modes  shown  in  Figure  l#x  is  negative# 
indicating  IniUblllty. 

There  are  several  interesting  features  of  Figure  1.  First# 
it  is  evident  that  the  growth  rates  of  the  modes  with  thermal 
boundary  condition#  Ti,ba88  *0#  are  insensitive  to  the  value  of 
the  base  temperature  or  to  the  mass  of  chromospheric  material. 
In  particular#  the  first  symmetric  mode  hasX  -  -.65;  the  first 


antisymmetric  mode  can  be  shown  analytically  to  have  1  ■  0  (Craig 
and  McClymont  1981).  The  modes  with  a  vanishing-perturbation 
heat  flux,  on  the  other  hand*  are  sensitive  to  the  value  of  the 
base  temperature#  in  agreement  with  our  previous  results 
(Antiouhos  1979);  but  only  if  the  amount  of  chromospheric 
material  is  very  small#  -  .1%  of  the  total  loop  mass.  For  models 
wxcn  a  significant  amount  of  cool  material#  the  growth  rates  are 
independent  of  the  thermal  boundary  conditions  or  of  the  amount 
of  base  material. 

We  conclude#  therefore#  that  the  static  models  of  coronal 
loops  are  linearly  unstable#  at  least  for  the  class  of  models 
considered  here.  The  fastest  growing  mode  in  all  cases  is  the 
lowest  order  one;  hence#  the  instability  is  a  global  one 
involving  the  whole  coronal  loop.  This  can  be  seen  in  Figure  2# 
where  we  have  plotted  the  eigenfunction  Tj(x)  for  the  case  with 
the  largest  chromosphere#  ~  20%  of  the  total  loop  mass.  Only  the 
mode  with  Tx.base  *  0  i8  ■hown*  however#  the  mode  with  T^base  " 
0  is  indistinguishable  from  it  except  very  near  the  base.  The 
modes  exhibit  a  significant  amplitude  only  in  the  coronal  and 
transition  region  part  of  the  loop.  This  amplitude  peaks  at 
approximately  the  temperature  where  the  radiative  loss  rate  n2  A# 
is  maximum  **  8  x  10*  K.  Below  this  temperature  the  amplitude 
falls  off  exponentially.  Although  the  peak  occurs  in  the  transi¬ 
tion  region#  the  coronal  contribution  to  the  "energy”  in  the 
mode#  i.e.#  /  Tj2  dx#  is  approximately  equal  to  the  transition 
region  contribution.  This  emphasises  the  point  that  the 
instability  is  a  global  one  and  not  restricted  solely  to  the 


transition  region 


B.  Rigid-Wall  Kodas 

In  Figure  1  X  is  plotted  Cor  the  first  syaaetric  rigid-wall 
■odes  and  Cor  each  of  the  two  theraal  boundary  conditions.  (The 
antisyaaetric  nodes  are  identical  to  those  of  the  free-surface 
case  as  discussed  above.)  It  is  evident  that  the  values  of  X  are 
very  different  froa  the  free-surface  case.  First,  they  are  all 
positive,  indicating  stability,  except  for  the  cases  with  very 
lictle  chromosphere,  <  .1%,  and  T^'  ■  0.  For  these  cases  the 
value  of  X  is  alnost  identical  to  the  corresponding  free-surface 
aoces.  However,  for  the  physically  relevant  cases  with  soae 
chroaospheric  aass,  the  values  for  X  are  positive.  Note  also 
that  for  nodels  with  significant  aaounts  of  chroaospheric  aass,  > 
10%,  the  values  of  X  decrease,  indicating  a  tendency  toward 
instability  as  the  chroaospheric  aass  increases. 

The  situation  can  be  clarified  by  inspection  of  the 
eigenfunction  Tj  (x)  for  the  lowest  syaaetric  aode  (Figure  2). 
As  in  the  free-surface  case,  the  aode  peaks  at  T  "  8  i  104  K,  but 
note  that  it  has  a  sero  crossing,  whereas  the  lowest  free-surface 
aode  can  have  none.  The  aaplitude  in  the  cbroaosphere  is  large? 
it  doainates  the  contribution  to  the  "energy”  in  the  aode.  As 
the  chroaospheric  aass  increases,  this  aaplitude  increases  and 
becoaes  aore  constant.  Bence  the  syaaetric  rigid-wall  aodea  are 
priaarily  chroaospheric  perturbations  rather  than  coronal  pertur¬ 
bations  as  in  the  free-surface  case.  Since  our  aodel  for  the 
chroaosphere  is  highly  artificial,  the  physical  significance  of 
the  rigid-wall  aodes  is  questionable. 
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The  key  result  of  this  paper  is  that  the  static  aodels  of 
coronal  loops  are  theraally  unstable,  even  if  the  aodels  contain 
arbitrary  aaounts  of  chroaospheric  aaterial.  This  disagrees  with 
the  conclusions  of  several  authors  (Habbal  and  Rosner  1979, 
Chiuderi,  Einaudi  and  Torricelli-Ciaaponi  1981,  Craig  and 
McClyaont  1981  and  McClyaont  and  Craig  1981a, b,c),  but  agrees 
witn  tne  conclusions  of  Hood  and  Priest  (1980)  and  Wragg  and 
Priest  (1982).  The  growth  rate  of  the  instability  is  of  the 
order  of  the  cooling  tiae  of  the  coronal  plasaa,  independent  of 
the  particular  paraaeter  of  the  aodel.  We  can  identify  two 
possible  physical  aanifestations  of  this  instability.  One  is  in 
the  observations  of  velocities  in  the  transition  region  (e.g., 
Lices  et  al.  1976).  Since  the  aaplitude  of  the  unstable  aode 
peaks  in  the  transition  region,  one  night  expect  the  largest 
velocities  to  be  generated  there.  Another  possible  nanifesta- 
tion  of  thernal  instability  is  in  observation  of  cool 
condensations  in  the  corona)  particularly,  quiescent  proainences. 
Our  results  inply  that  coronal  loops  are  naturally  unstable)  and, 
hence,  condensations  aay  fora  spontaneously  in  the  corona. 

Of  course,  in  order  to  investigate  these  possibilities,  the 
nonlinear  evolution  of  the  instability  aust  be  calculated. 
Clearly,  the  instability  aust  not  saturate  at  a  low  level  if  it 
is  to  produce  sisable  velocities  in  the  transition  region  or  cool 
B-alpha  condensations  in  the  corona.  As  discussed  in  our 
previous  paper  (Antiochos  1979),  we  believe  that  the  instability 
will  result  in  a  nonlinear  oscillation  of  the  loop  about  its 
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static  stats.  Ths  period  of  tbs  oscillation  should  be  related  to 
the  growth  rate  of  the  instability,  103  sec.  for  typical  loop 
parameters.  The  amplitude  of  the  oscillation,  which  must  be 
calculated  by  nonlinear  simulation,  will  determine  the  physical 
importance  of  the  instability. 

Numerical  simulations  of  coronal  loops  have  been  performed 
by  several  authors.  These  simulations  have  not  exhibited  any 
significant  evidence  for  thermal  instability.  However,  in  all 
cases  rigid-wall  boundary  conditions  were  used,  and  it  is  clear 
from  our  results  here  that  these  conditions  will  tend  to 
stabilize  the  models,  at  least  linearly.  In  addition,  numerical 
effeccs  may  be  providing  unphysical  stability  to  the  models.  We 
intend  to  consider,  in  detail,  the  nonlinear  stability  of  the 
static  loop  models  in  a  forthcoming  paper. 
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flams  Captions 

Ei8U££  1 

Absolute  value  of  *  for  a  series  of  static  aodels  with 
decreasing  base  teaperature  and/or  increasing  chronospheric  aass. 
The  results  for  five  aodes  are  shown.  Solid  lines  refer  to  free* 
surface  aodes,  dashed  lines  refer  to  rigid-wall  nodes. 

flame  2 

The  eigenfunction  Tj  for  the  lowest  synaetric  free-surface 

(solid  line)  and  rigid-wall  (dashed  line)  aodes.  They  are 

plotted  as  a  function  of  log  (r)  where  r  is  defined  ass  r  • 
x 

/n~(x)  x  dx  . 


